Shot noise and multiple Andreev reflections in d-wave superconductors 
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We present a theoretical analysis of the shot noise in d-wave/d-wave contacts with arbitrary 
transparency, including the contribution of multiple Andreev reflections. The multiple charge quanta 
transferred in these processes are revealed as a huge enhancement of the noise-current ratio at low 
voltages, which survives for all crystal misorientations. We also show how different ingredients like 
non-magnetic impurities or a magnetic field produce very characteristic hallmarks in the shot noise, 
which can be used as a further test of the d-wave scenario in superconducting cuprates. 
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In the last years the intensive study of the nonequilib- 
rium current fluctuations, known as shot noise, has pro- 
vided a deeper understanding of the electronic transport 
in many different contexts . The shot noise reveals as- 
pects hidden in the usual conductance measurements like 
the statistics and charge of the carriers, relevant energy 
scales or transmission information. In the case of su- 
perconducting contacts, the noise has been mainly used 
for the analysis of the effective charges transferred in the 
different tunneling processes. This effective charge can 
be roughly defined as the noise-current ratio. A strik- 
ing example is the recent observation in superconducting 
atomic-size contacts of effective charges much larger than 
unity attributed to multiple Andreev reflections S, and 
in quantitative agreement with the theoretical expecta- 
tions ||. Unfortunately, the analysis of shot noise has 
been mainly restricted to conventional superconductors, 
and only a few theoretical works have recently addressed 
noise in NIS-junctions with d-wave superconductors Q. 

On the other hand, the origin and nature of the high 
temperature superconductivity in the cuprates is still 
an open problem. Different phase-sensitive experiments 
have provided strong indications that the order parame- 
ter in these materials has a dominant d x 2_ y i component 
H Sj. However, these experiments have not definitively 
closed the debate about basic questions like the univer- 
sality of this symmetry or the existence of subdominant 
components [0, §]. In this sense, it is highly desirable 
to provide new tools which can further test the differ- 
ent scenarios. Thus, a natural question is: what can the 
shot noise teach us in unconventional superconductivity? 
For instance, the shot noise provides fundamental infor- 
mation on the charge of the carriers, as it was shown in 
complex situations like in the fractional quantum Hall 
effect Q or in superconducting point contacts [0. On 
the other hand, the shot noise in a junction depends in 
a different way on the interface properties as compared 
with the current. In this sense, the study of this quan- 
tity can used as a cross-check for the different transport 
theories, and in turn it can be very valuable to solve the 
lack of consensus in the interpretation of the tunneling 



experiments in cuprate junctions ]10[ . 

In this paper, we present the first theoretical analy- 
sis of the shot noise in d-wave/d-wave SIS junctions of 
arbitrary transparency. We shall show that the zero- 
frequency noise, S, may by large exceed the Poisson value 
2el, where / is the current, due to the occurrence of mul- 
tiple Andreev reflections. In particular, at high trans- 
parencies the effective charge, q, defined as q — S/2I, 
exhibits a huge enhancement at low voltages (q/e 3> 1), 
which survives for all crystal misorientations. At low 
transparencies, contrary to the s-wave case, q is not quan- 
tized in units of the electron charge due to the averaging 
over the anisotropic gap. We shall also show that elastic 
scattering mechanisms like bulk impurities may result in 
a strong reduction of the effective charge. Finally, we 
shall show how the Doppler shift of the Andreev bound 
states in the presence of a magnetic field is revealed in 
the shot noise. All these features are very characteristic 
of the d-wave symmetry and can be used as additional 
tests of this scenario in cuprates. 

Our goal is to extend the theory of the shot noise to 
the case of superconducting cuprates. For this purpose, 
we consider a voltage biased contact, consisting of two 
d x 2_ y 2 superconductors separated by a single interface 
of arbitrary transparency. The order parameter on side 
i, i = L, R, is rotated by on with respect to the sur- 
face normal, and we denote junction type by the rela- 
tive crystal orientations as d aL -d aR . There are several 
experimental realizations of this system, among which 
the bicrystal grain-boundary junctions are ideal exam- 
ples |pj| . To calculate the noise we use the formalism 
developed in Ref. |Q. In that work we introduced a 
formulation of boundary conditions that mimics inter- 
faces for the quasiclassical theory of superconductivity 
and that are suitable for arbitrary transparency, and 
we established the machinery to determine the current 
fluctuations in unconventional junctions. Here we con- 
sider the case of point-contact-like geometry and assume 
that the voltage drop takes place at the interface. Thus, 
to compute the noise we first determine self-consistently 
the local electronic properties of the isolated electrodes. 



2 



This includes effects on the order parameter profile and 
on the local density of states (DOS) by pair breaking 
caused both by quasiparticle scattering off the interface 
and off homogeneously distributed impurities in the crys- 
tals jl3], Q 0, Finally the noise is calculated 
using the local surface Green's and solving the appropri- 
ate boundary conditions for a point contact, as detailed 
in Ref. Q. 

The noise spectral density S(lu) is defined as 



S(u) = / d(t' - t) e tuJ ^'- t ){5i{t l )5i{t) + 6i(t)SI(t')), 



where 5I(t) = I(t) — (I(t)) are the fluctuations in the 
current. We only consider the zero- frequency limit at 
zero temperature. In the case of a constant bias volt- 
age, V, one can show (see Ref. [|l2|) that the noise os- 
cillates in time with all the harmonics of the Joseph- 
son frequency, i.e. S(t) = J2m SmC 1 ™ 1 ^, where <f>(t) — 
4>a + (2eV/h)t is the time-dependent superconducting 
phase difference. We shall only consider the dc noise, de- 
noted from now on as S. Furthermore, we assume that 
the interface conserves the momentum of the quasiclas- 
sical trajectories, which allows us to write the noise as 
a sum over independent trajectory contributions: S — 
\ /-tt/2 ^ S(i^)cos(%), where defines the Fermi 
surface position. For the angular dependence of the 
transmission coefficient we use the expression D(q.) — 
D cos 2 (££,)/ [1 — Dsai (fy.)], resulting from a (5-like poten- 
tial. Here D is the transmission for the trajectory perpen- 
dicular to the interface. In the tunneling regime one can 
easily demonstrate that the zero-frequency noise reaches 
the Poisson value, i.e. S — 2el. Thus, in this limit 
the noise does not contain new information as compared 
with the current. For this reason, we shall investigate 
the case of not too low interface transparency, D > 0.1, 
in which the multiple Andreev reflections (MAR) play a 
fundamental role in the transport []l9[ pp[ . 

Let us start by analyzing the case of a symmetric do- 
do junction in the clean limit. In this case, the order 
parameter is constant up to the surface, and there are no 
bound states for any trajectory. The noise-voltage char- 
acteristics for a single trajectory, S(ffc), coincide with 
those of isotropic s-wave superconductors || , and can be 
seen in Fig |l|(a). As a consequence of the occurrence 
of MARs, the trajectory-resolved shot noise exhibits the 
following remarkable features: (i) the presence of a pro- 
nounced subharmonic gap structure (SGS) at voltages 
eV = 2A(^.)/rt, (ii) the noise greatly exceeds the Pois- 
son value 2el in the subgap region, as can be seen in 
Fig |l](b), and (iii) in the tunneling regime the effec- 
tive charge is quantized in units of the electron charge. 
This last feature, illustrated in the inset of Fig |l|(b), was 
used to suggest that the noise provides a way of mea- 
suring the charge of individual MARs in s-wave super- 



conductors ||. Indeed, as mentioned in the introduc- 
tion, these noise-voltage characteristics have been quan- 
titatively confirmed in the context of superconducting 
atomic contacts 0|. The natural questions now are: do 
these features survive after doing the average over the 
different directions in the Fermi surface? Can we still 
identify the charge of individual MARs in a d-wave junc- 
tion?. The answers to these questions can be seen in 
Fig. [l](c-d). First of all, notice that the SGS is still visi- 
ble, but it is more rounded than in the s-wave case. It is 
worth remarking that it is the bulk maximum gap what is 
revealed in the SGS. Notice also that the effective charge 
does not show any sign of quantization even at low trans- 
parencies (see inset of Fig. 0(d)). This is due to the fact 
that different MARs contribute simultaneously for differ- 
ent trajectories, and then the discreteness of q is washed 
out. Anyway, the dominant contribution of MAR at high 
transmission is still manifested as a huge enhancement of 
the effective charge at low bias (q e). 
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FIG. 1: do-do contact in the clean case: (a) Angle-resolved 
shot noise as a function of the voltage in units of the trajec- 
tory gap for different transmissions. Gjv is the normal state 
conductance, (b) Angle-resolved effective charge, q = 5/27. 
(c) Angle-averaged shot noise. The voltage is normalized by 
the maximum gap A. (d) Angle- averaged effective charge. 



Let us now consider the case of a d 7T /4~d_ 7T /4 junction. 
In this case, assuming specular quasiparticle scattering 
at the interface, an Andreev bound state forms at zero 
energy for every trajectory pi[ . This implies that the 
surface acts as a pair- breaker ||, and the gap is 
depressed in the vicinity of the interface, vanishing ex- 
actly at the barrier. This order-parameter profile induces 
not only the appearance of bound states at zero energy, 
but also at the gap edges, as can be seen in the inset of 
Fig. ||(b). As a consequence of this local density of states 
the noise exhibits a pronounced SGS due to resonant tun- 
neling through the bound states (see Fig. 0(a)). As in the 
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FIG. 2: 
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_ x /4 contact in the clean case: (a) Shot noise 
as a function of voltage for different transmissions, (b) Ef- 
fective charge vs voltage. The curves were computed using a 
small inelastic broadening (~ 0.003A). Inset: local DOS at 
the interface for a 45° misorientation for different values of 
the bulk-impurity scattering rate F (Born scatterers), mea- 
sured in units of 2nTc, where Tc is the critical temperature 
in the clean case. A is the maximum bulk gap for the clean 
superconductor. 



case of the current, see Ref. |2Cj, there is an even-odd ef- 
fect in the SGS, in the sense that the even (eV = A/n) 
structures are more pronounced. Its origin can be under- 
stood as follows. In this geometry there are two types 
of MARs which dominate the transport: (a) those which 
connect the bound states with the gap edges and (b) the 
usual ones connecting the gap edges. The first ones give 
rise to the SGS at eV — A/n, while the second ones 
contribute to the whole series eV — 2 A/n. However, the 
bound states at the gap edges do not appear for every tra- 
jectory, which weakens the SGS due to these processes. 
On the other hand, as we show in Fig. ||(b) the effec- 
tive charge is not even quantized at low transparencies, 
again due to the average over the different trajectories. 
However, at high transparencies the dominant contribu- 
tion of the MARs give rise to a huge enhancement of 
the effective charge at low bias. This is a robust feature 
which survives for all crystal misorientations, and it is an 
unambiguous signature of the fact that the MARs con- 
trol the low voltage transport. Indeed, this pronounced 
increase of the noise-current ratio at low bias has been re- 
cently observed in symmetric bicrystal YBCO Josephson 
junction [p2| in, to our knowledge, the first experimental 



analysis of the shot noise in cuprate Josephson junctions. 
In this experiment a mean transparency of D « 0.01 was 
estimated, but in our opinion this enhancement is due to 
MARs in high transparent conduction channels, probably 
due to the presence of pinholes like in the conventional 
SIS tunnel junctions of Dieleman et al. p3[ . 

In d-wave superconductors the order parameter is very 
sensitive to scattering from nonmagnetic impurities and 
surface roughness. In particular, it is known that these 
elastic scattering mechanisms provide an intrinsic broad- 
ening for the zero-energy bound states [jl7|, |igfl. F or the 
case of Born scatterers this broadening is cx vTA, where 
r = l/2r is the effective pair-breaking parameter locally 
at the surface. This is illustrated in the inset of Fig. ||(b) 
for the case of bulk impurities. The interesting ques- 
tion now is: what is the signature of impurities in the 
shot noise of a d-wave junction? In Fig. || we show the 
shot noise and effective charge for a d^u-d^u junction 
for different values of the bulk-impurity scattering rate. 
As the scattering rate increases, there are two major ef- 
fects that one should notice: (i) the disappearance of the 
SGS in the noise, and (ii) a reduction of the effective 
charge, specially pronounced at low voltages. Both fea- 
tures can be understood as follows: the increase of den- 
sity of states in the gap region enhances the probability of 
single-quasiparticle processes, producing the subsequent 
reduction of the probability of the Andreev processes, 
which in turn leads to both the suppression of the SGS 
and the reduction of the effective charge. 
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FIG. 3: Shot noise and effective charge as a function of voltage 
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for a d^/4-d- 7r /4 contact for two values of the bulk-impurity 



Fogelstrom et al. [[L5[ have shown that the Andreev 
bound states should split in the presence of a magnetic 
field perpendicular to the a6-plane. This splitting results 
in a splitting of the zero bias conductance peak observed 
in tunnel junctions j24j. It is then interesting to ana- 
lyze what is the signature of this time reversal symmetry 
breaking in the shot noise. Let us consider a magnetic 
field perpendicular to the afc-plane, H = Hz. As men- 
tioned above, this field leads to a Doppler shift in the 
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FIG. 4: Differential shot noise of a clean d^i^-d-^/A junction 
with D = 0.2 and different values of the magnetic field. The 
curves has been vertically displaced for clarity and dotted 
lines have been added to guide the eye. 



continuum excitations given by ■ p s , where the con- 
densate momentum is p s = —{e/c)A(x)y, with A the 
self-consistently determined vector potential jl5|. This 
means that the Andreev bound states are shifted to an 
energy which, in the limit of a large ratio A/£o, can be 
estimated to be eb{%) = {e/c)VfHXsxsxj^, A being the 
ab-plane penetration depth. We shall use a natural field 
scale set by a screening current of order the bulk critical 
current, Hq — cA/eu/A, which is of the order of a Tesla 
p5[ . The screening currents flow parallel to the inter- 
face and in opposite directions in both electrodes, which 
means that the trajectory resolved DOS of the left and 
right superconductors are shifted by 2t b {%) relative to 
each other. As explained in Ref. |^(|, this shift modifies 
the threshold voltages of MARs starting and ending in 
different electrodes, leading to the splitting of the peaks 
with an odd order n in the SGS. On the contrary, since 
the magnetic field produces a rigid shift of the spectrum, 
the threshold voltages of those MARs starting and end- 
ing in the same electrode are not modified. This means 
that the positions of the structures with an even order 
n in the SGS remain unchanged. This is illustrated in 
Fig. |^ where we show the differential shot noise, dS/dV, 
for a d n u-d_ w /4 junction with transmission D = 0.2 for 
different values of the magnetic field. Starting at large 
voltages, the structure at 2A splits with applied field. 
Around eV = A there is a maximum at eV = A, un- 
affected by the applied field, as well as a field-shift of 
the peak just above A. The field dependence of the dif- 
ferential noise is most clearly resolved at larger biases, 
eV > A/2, as the marks of the various processes begin 



sin ~ and thus having a vanishing Doppler shift. 

In summary, we have presented a theory of shot noise 
in d-wave/d-wave contacts with arbitrary transparency. 
We have shown that in the MAR regime these nonequilib- 
rium current fluctuations exhibit very peculiar features 
like subharmonic gap structure, super-Poissonian noise 
(S 3> 2el), reduction of the effective charge q by impu- 
rities and the splitting of the SGS in a magnetic field. 
All these features are unique fingerprints of the d-wave 
scenario, and we hope that our analysis will trigger off 
experimental study of shot noise in cuprate junctions. 
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Dynamics, and MF by the Swedish research council. 
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